Collisional stability of a three-component degenerate Fermi gas 
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We report on the creation of a degenerate Fermi gas consisting of a balanced mixture of atoms in 
three different hyperfine states of 6 Li. This new system consists of three distinguishable Fermions 
with different and tunable interparticle scattering lengths 0,12, 0,13 and <J23- We are able to prepare 
samples containing 5 ■ 10 4 atoms in each state at a temperature of about 215 nK, which corresponds 
to T/Tf ~ 0.37. We investigated the collisional stability of the gas for magnetic fields between 
and 600 G and found a prominent loss feature at 130 G. From lifetime measurements we determined 
three-body loss coefficients, which vary over nearly three orders of magnitude. 

PACS numbers: 



Since the first successful preparation of a degenerate 
Fermi gas of 40 K atoms in 1999 [l|, tremendous progress 
has been made in the field. The key to this success 
was the possibility to tune the interactions in mixtures 
containing two spin components at wish using Feshbach 
resonances Large, resonant scattering lengths facili- 
tated evaporative cooling and led to the first observation 
of a strongly interacting Fermi gas Q. At large, posi- 
tive scattering lengths, weakly bound molecules could be 
formed [J] . These bosonic molecules were condensed into 
a molecular Bose-Einstein condensate (BEC) [1, @, 0]- 
The tunability of the scattering length across the Fesh- 
bach resonance gave experimentalists access to the so- 
called BEC-BCS crossover , which connects in a 
continuous, isentropic way molecular Bose-Einstein con- 
densates with two-component, BCS-type Fermi gases. In 
this limit the attractive interaction causes the atoms to 
form weakly bound many-body pairs pdj . In the inter- 
mediate range, where the scattering length diverges, the 
so-called unitary regime could be accessed [T3|, where 
the interactions depend only on the inverse Fermi mo- 
mentum 1/kF, which causes the physics to become inde- 
pendent of any length scale and therefore universal for 
any Fermi system with infinite scattering length. Conse- 
quently, such tunable two-component ultracold systems 
can be an ideal model for such diverse systems as neutron 
stars [lH or, with the addition of an optical lattice [3], 
high-Tc-superconductors [l5[ . 

When a third spin state is added to such a two- 
component gas, a number of intriguing questions arises 
immediately, which have already been extensively stud- 
ied in theory [l6|, Qj], ■ How will pairing occur in such 
a system: Will the individual components compete, and 
only two of them form pairs, while the third component 
remains a spectator, or will the lowest-energy state of the 
system be a three-body bound state [H, US HH ? There 
are predictions for a phase transition between a super- 
fluid and trionic phase in optical lattices, which can be 
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treated analogous to baryon formation in QCD [22] . Fur- 
thermore, as the Fermi pressure is lower than in a two- 
component mixture of equal density, the stability of such 
a gas in the case of resonant two-particle interactions is 
still a controversial topic [23|, [lij . Ultracold ensembles of 
6 Li offer the unique possibility to study these phenomena, 
as wide, overlapping Feshbach resonances between the 
three lowest-lying hyperfine states allow to tune the two- 
body scattering lengths over a wide range (see Fig.[2]c)). 

So far there has been only little experimental work on 
three-component Fermi gases: A third non-degenerate 
component was used as a probe for thermometry [25| . In 
other experiments radio frequency (RF) transitions to a 
third state were used to study pairing in a two-component 
system (TTI . [26j, but a stable three-component mixture 
has not been observed. Here we present the first real- 
ization of a stable and balanced three-component Fermi 
gas, which is a first step towards exploring the rich phase 
diagram of this system. Furthermore, we investigate the 
magnetic field dependence of the collisional stability of 
this gas. 



900 






6> 






5> 








4 


600 








300 


- F=3/2 ^-^^^ 















300 




^^^^^^ applied radio 




600 














3> 








2> 


900 






1> 



m = 3/2 
m = 1/2 
m = -l/2 



m = -3/2 
5= -1/2 
m„= 1/2 



100 200 300 400 500 600 

magnetic field [G] 



FIG. 1: Zeeman hyperfine levels of Li in the electronic 
ground state. Transitions between adjacent hyperfine states 
can be driven with external RF fields. 

Our approach to the preparation of three-component 
Fermi mixtures makes use of well-established techniques 
used for studying ensembles of ultracold atoms [27j • We 
first precool 6 Li atoms in a magneto-optical trap (MOT) 
and transfer about 10 6 atoms into a far red-detuned op- 
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tical dipole trap. During the transfer the atoms are op- 
tically pumped into the F = 1/2 hyperfine state so that 
at this stage of the experiment only the two lowest mag- 
netic subleves (labeled |1) and |2) in Fig.[T|) are occupied. 
The dipole trap is formed by crossing two counterpropa- 
gating beams with a waist of 50 /mi, which results in a 
cigar shaped trap with an aspect ratio of about 10 : 1. 

The wide Feshbach resonance in 6 Li allows us to tune 
the two-body s-wave scattering length an between states 
|1) and 1 2) to 3560 ao at a magnetic field of 751 G, which 
leads to high thermalization rates. This enables us to 
perform fast and efficient evaporative cooling by lower- 
ing the laser power in the trap. The large positive value 
of the scattering length is associated with a bound state 
[28| , whose binding energy depends on the magnetic field. 
At 751 G the binding energy is fee • 2.3/iK, where ks is 
Boltzmann's constant. As soon as the temperature of 
the sample becomes comparable to this binding energy 
the bound state is populated through three-body recom- 
bination. By further evaporation these molecules can be 
condensed into a Bose-Einstein Condensate containing 
1.5 • 10 5 molecules. 

The third spin state (|3)) can be populated by driv- 
ing RF transitions between the different magnetic sub- 
levels, which has been extensively used to perform RF 
spectroscopy on two-component Fermi gases. In those 
experiments, which were performed close to the Fesh- 
bach resonances, high inelastic loss prevented the prepa- 
ration of stable three-component mixtures in 6 Li [26j. 
We overcome the problem of loss during preparation of 
the mixture by tuning the magnetic field close to the 
zero-crossings of the two-particle scattering lengths. To 
avoid populating the molecular state during evaporation, 
we tune the scattering length to a small negative value 
of a\2 ~ — 300 ao at B = 300 G before molecule forma- 
tion sets in. There we continue evaporative cooling of 
the two-component mixture, reaching a temperature of 
130 nK (T/T F w 0.28). We then minimize the interac- 
tions by ramping the magnetic field to a value of 563 G. 

In order to prepare an equal mixture of atoms in states 
|1), 1 2) and [3} we simultaneously drive the RF transi- 
tions |1) <-> \2) and \2) <-» |3) for a period of 850ms at a 
magnetic field of 563 G. A combination of collisions and 
a residual magnetic field gradient lead to an incoherent 
mixture of atoms in the three lowest spin states. We can 
produce samples containing 5 TO 4 atoms per spin state at 
a temperature of 215 nK in a harmonic trap with trap fre- 
quencies uj x — ujy = 27r-386(15) Hz and oj z — 27r-38(2) Hz. 
Therefore, the starting point for our experiments is a de- 
generate three-component Fermi gas with T/Tp rs 0.37 
and a peak density of 6 • 10 11 atoms/cm 3 . Under these 
conditions the gas is stable with a 1/e lifetime greater 
than 30 seconds. 

In our first experiment we studied the magnetic field 
dependence of the collisional stability. We prepared a 
three-component mixture and held the atoms at vari- 
ous magnetic fields between and 750 G for 250 ms. To 
record the remaining fraction of atoms in one of the spin 



states selectively, we make use of their Zeeman splitting. 
We calibrated the imaging laser frequencies to the respec- 
tive transitions at a magnetic field of 526 G. Therefore we 
jump to this field at the end of each experimental cycle 
and record the number of atoms in one spin state via ab- 
sorption imaging. The result is shown in Fig. [2] a). Each 
data point is the average of 5 individual measurements. 
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FIG. 2: a) Magnetic field dependence of the remaining frac- 
tion of atoms after holding a three-state mixture for 250 ms. 
b) Result for a two-state mixture under the same experi- 
mental conditions. The enhanced trap loss at fields larger 
than 600 G is due to high inelastic collision rates close to the 
two-body Feshbach resonances, which also causes loss in the 
three-component mixture. For magnetic fields larger 600 G, 
molecules are formed in the |1) - |3) channel which are not 
detected when ramping back for imaging, c) Two-body scat- 
tering lengths for all spin-state combinations 29] (an = for 
symmetry reasons). 

In the region from 540 G to 590 G, where all 
two-particle scattering lengths cross zero, the three- 
component mixture is stable and the remaining fraction 
is close to one. For larger magnetic fields one observes 
a strong inelastic decay of the sample due to the two- 
body Feshbach resonances. Below 540 G the fraction of 
remaining atoms quickly drops to a minimum of 0.6 at 
470 G and then increases linearly to 0.75 at 200 G. Be- 
tween 80 G and 190 G there is a strong loss feature where 
the fraction of remaining atoms drops to about 0.5 at 
130 G. For magnetic fields below 70 G, where ai2 and 0,23 
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are smaller than 60 do, the mixture is once again stable. 

In order to separate processes involving only two dis- 
tinguishable Fermions, we did reference measurements 
for all possible two-state mixtures by preparing three- 
component mixtures and removing one component with 
a resonant laser pulse before ramping to the probe field. 
As an example, the result for the mixture of atoms in 
states |1) and |3) is shown in Fig.[5]b). Away from the 
Feshbach resonances all two-component mixtures are sta- 
ble due to Pauli blocking, as a three-body loss process 
would involve at least two identical Fermions. Adding 
a third distinguishable Fcrmion allows three-body pro- 
cesses, which can decrease the stability of the gas. From 
this we infer that all decay observed in Fig. [2] a) for mag- 
netic fields below 590 G stems from processes involving 
atoms in all three spin states. 

To get quantitative information about the loss rates 
we measured the decay of the gas over a period of five 
seconds at several magnetic field values and hence for 
different two-body scattering lengths. For these mea- 
surements we prepared the three-component mixture as 
described above, tuned the magnetic field to the value 
of interest and measured the number and temperature of 
remaining atoms as a function of time. As an example, 
Fig. [3] shows the decay curve for state |2) at a magnetic 
field of 300 G. We observed that the ratio of particles 
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FIG. 3: Evolution of a) atom number and b) temperature of 
atoms in state |2) over five seconds at 300 G. T denotes an 
effective temperature deduced from a gaussian fit to the cloud 
after time of flight. Each data point is a mean value of three 
independent measurements taken in random order. The solid 
line is a fit to the data, applying the method described in the 
text. 

in the three states remains one during the whole decay. 
Therefore we can assume in the following analysis that 
all states have equal loss rates, that are governed by a 
three-body process. The loss of particles can then be 
described by 

h l (r) = -K 3 n l (r) 3 , (1) 

where is the three-body loss coefficient and rii{r) de- 
notes the local, temperature dependent density of atoms 
per spin state \i). Additionally, we take into account a 
small one-body loss with a 1/e-lifetime of about 100 s. 
From the two-component measurements we infer that 
we can neglect two-body decay processes. To deduce 
K3 from our experimental data, we follow a method de- 
scribed in [30l | . Loss of atoms occurs predominantly in 



the center of the trap where the density is high. The en- 
ergy taken away by a lost particle is thus smaller than the 
mean energy per particle when averaging over the whole 
cloud. This leads to heating of the sample (see. Fig.[3]b) 
which requires a numerical treatment of the atom num- 
ber and temperature evolution, for which we use the code 
developed for [3(|. The model we use describes loss in 
a thermal gas. For T/Tp ~ 0.37 minor changes occur 
due to degeneracy, which we neglect for simplicity. For 
the analysis of the lifetime curves we use an effective 
temperature (T in Fig. [3]) deduced from a gaussian fit 
to the density distribution after time of flight. If the 
gas is degenerate, the obtained value T is slightly higher 
than the real temperature. As the density in a degener- 
ate Fermi gas is reduced with respect to a thermal gas 
of the same temperature, the higher temperature value 
compensates to some extent for the small effects of de- 
generacy. Hence, we approximate T k,T for our analysis. 
For most values of K 3 this approximation affects only the 
very first data points, as the temperature quickly exceeds 
the Fermi temperature after initial loss and heating. The 
solid line in Fig. [3] shows the fit to the atom number and 
temperature evolution according to this model for one 
decay curve. From such fits we obtain our values of K3. 

Fig.|4] shows the obtained three-body loss coefficients 
for all three species as a function of the magnetic field. 
As expected, K3 reflects the qualitative behavior of the 
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FIG. 4: Three-body loss coefficient K3 vs. magnetic field. 

atom number in Fig. [2] The relative error of the data 
points is caused by the uncertainty in spin balance, dif- 
ference in detection efficiency for the three states and 
fitting uncertainties. This error can be estimated by the 
scatter of the fitted values for K3 for different states and 
is much smaller than the observed variation of K 3 . The 
absolute scale is subject to additional systematic errors 
in particle number and trap frequency. As K3 oc N 2 u) 6 
uncertainties of 40% in N and 7% in uj lead to an error 
of 90% in the absolute scale. 

Although there is currently no theoretical model which 
can quantitatively explain the strong variation of K3 with 
the magnetic field, we can qualitatively discuss some as- 
pects of its behavior. In a naive picture one can describe 
a three-body event by elementary two-particle interac- 
tions governed by ai2, 013 and 023. If two of the three 
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scattering lengths are close to zero, only two of the three 
spin states are interacting while one remains a spectator. 
Then the physics is comparable to a two-component sys- 
tem and the gas is stable. If only one scattering length 
dij is small, interaction between states \i) and \ j) can be 
mediated by the third particle, allowing three-body pro- 
cesses. This is consistent with our observation that in the 
regions below 70 G, where a\2 and 023 are close to zero 
and around 570 G, close to the zero-crossings of a 13 and 
023, -K3 becomes minimal. The smallest observed value 
of 10~ 25 cm 6 /s is still four orders of magnitude larger 
than the value for three-body recombination in a thermal 
cloud of 87 Rb [31|, but still sufficiently small to observe 
lifetimes on the order of 30 s at our densities. 

At the loss feature near 130 G the values of K3 are 
enhanced by three orders of magnitude compared to the 
minimal value. In this magnetic field range all two-body 
scattering lengths are negative, consequently no weakly 
bound dimer state exists in this field range. As all three 
scattering lengths are large compared to the effective 
range of the interatomic potential, one would expect its 
short-range properties not to be relevant and therefore 
the physics to be universal. If this were the case the ob- 
vious explanation for the loss feature would be a trimer 
state crossing into the continuum, although we cannot 
exclude other possible scenarios such as relaxation into 
deeply bound dimer states or spin changing collisions. 
The properties of such trimers would be comparable to 
Efimov trimers [32] , which were observed through an en- 
hanced three-body loss rate in bosonic Cs atoms [33j ] - 
An essential difference is that our system consists of dis- 
tinguishable Fermions, therefore the formation of such 
trimers would have almost the same symmetry as the 



formation of baryons from three quarks in QCD. Further 
theoretical and experimental investigation will be nec- 
essary to quantitatively understand the magnetic field 
dependence of the three-body decay rate. 

In this work we have established a method to create a 
stable degenerate Fermi gas containing atoms in three 
different spin states of 6 Li. We have determined the 
three-body decay rate of the gas for magnetic fields be- 
low 600 G, whose behavior we can in part explain with 
the two-body scattering lengths between the atoms. The 
most prominent observation is a broad loss feature caused 
by a process involving atoms in all three spin states, 
which indicates the presence of a three-particle reso- 
nance phenomenon. The preparation of a stable three- 
component gas in thermal equilibrium opens up new pos- 
sibilities to study many-body phenomena in ultracold 
gases. As an example, one could switch the scatter- 
ing lengths rapidly to large values and observe the con- 
sequent collective evolution of the cloud, which would 
reflect the equation of state of the many-body system. 
At sufficiently low density, the collective time evolution 
should occur at a faster time scale than the loss of atoms 
from the trap. 
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Recently we have learned that in the group of 
K.M. O'Hara similar measurements have been performed 
independently [34| . 
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